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We study the effects of the Coulomb interaction in the one dimensional Kondo lattice model on 
the phase diagram, the static magnetic susceptibility and electron spin relaxation. We show that 
onsite Coulomb interaction supports ferromagnetic order and nearest neighbor Coulomb interaction 
drives, depending on the electron filling, either a paramagnetic or ferromagnetic order. Furthermore 
we calculate electron quasiparticle life times, which can be related to electron spin relaxation and 
decoherence times, and explain their dependence on the strength of interactions and the electron 
filling in order to find the sweet spot of parameters where the relaxation time is maximized. We find 
that effective exchange processes between the electrons dominate the spin relaxation and decoherence 
rate. 
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I. INTRODUCTION 



Recently, the interest in nanoscale systems has been 
rapidly increasing. Among them are 13 C carbon 
nanotubeSfi^ nanowires^i and carbon nanotubes filled 
with cndohcdral fullercncs or molecular magnets^. The 
above mentioned systems have in common, that they 
consist of local spins (electron or nuclear spins) which 
interact via exchange interaction with itinerant conduc- 
tion electrons. These are exactly the constituents of 
the one dimensional Kondo lattice model&i (KLM). To 
make these materials available for spin electronics or 
quantum information processing it is necessary to un- 
derstand their properties in detail: ground state (e.g. 
magnetic order), spectral (e.g. dispersion relation of elec- 
trons) and dynamical (e.g. non-equilibrium, spin relax- 
ation/decoherence) properties. 

Interaction between the local spins in the KLM is 
generated effectively due to the hopping t of electrons 
and an onsite direct spin exchange J between the itin- 
erant and localized spins, see Fig. [TJ This interaction 
is a result of the competition of onsite singlet forma- 
tion and an effective RKKY (Ruderman-Kittel-Kasuya- 
Yosida) interaction^. The order of the local spins due to 
the interaction is captured in the phase diagram of the 
KLM,— which is basically divided into three phases 
depending on J/t and the electron filling n (n = 1 is half 
filling). At n = 1 the system turns out to order anti- 



ferromagnetically for arbitrary coupling strength. A fer- 
romagnetic (FM) phase is established, if either J is large 
enough or n is small enough^ Otherwise the local spin 
lattice is in the paramagnetic (PM) phase, because then 
the effective RKKY interaction dominates the system. 

The mechanism of ferromagnctism in the KLM can 
also be understood in terms of an electron quasiparticle 
picture, where the quasiparticle is the so called spinpo- 
laro n 15 ' 16 , see Fig. [2^. For a given FM order of the lo- 
cal spins in a Id system it was shown that the itinerant 
electrons and the magnons of the local spin bath form a 
bound spinpolaron state which is detectable in transport 
measurements and was proposed as a long-living corre- 
lated many-body spin stated formi ng p ossibly one part 
of a many-body spin qu-bit. In Ref. LLJ it was shown for 
the case of a single conduction electron that a spinpo- 
laron develops with a huge extent over the whole lattice 
leading to FM order in the ground state. In Ref. [13 this 
was extended to finite electron fillings and it was shown 
that long quasiparticle life times are connected with FM 
order of the local spins. In Ref. [3, the quasiparticle dy- 
namics of the half filled KLM (n=l) have been examined 
as well. By means of a strong coupling expansion up to 
11th order it has been possible to calculate the quasipar- 
ticle dispersion relation to good accuracy and it could be 
shown that the quasiparticlcs behave like nearly localized 
f-electrons due to the strong correlation of the conduction 
and localized electrons. 

It is known that the main relaxation and decoherence 
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FIG. 1. (Color online) The Kondo lattice model. The con- 
duction electrons are depicted in the upper row (red) and the 
localized electrons are depicted as bold arrows in the lower 
row (grey). 



from being numerically exact, acting in the lattice space 
without any approximations and taking all backaction 
effects of the local spin lattice on the conduction elec- 
trons automatically into account. Furthermore it allows 
for calculations in a broad parameter regime and works 
especially well for one dimensional systems with open 
boundary conditions and finite lattices. Here we are par- 
ticularly interested in finite lattices, since nanoscale sys- 
tems have finite sizes and show corresponding effects. 



source of single electron spins in semiconductor based 
quantum dots arises from interactions with the nuclear 
spin background i 18 i 25 ' 26 An appropriate path to diminish 
the relaxation is the application of a large magnetic field, 
whereas the decoherence rate is reduced by state distri- 
bution narrowing] 19 i 36 However, the initial preparation 
of the nuclear bath in a pure state (e.g. full polariza- 
tion) is an experimental challenge. Recently, the idea was 
proposed to consider the nuclear bath at very low tem- 
peratures in the FM phase, which is mediated by many 
itinerant electons via the RKKY-interaction^iiii 2 ^ In 
Ref.i it was shown that the Coulomb interaction in a 2d 
electron gas leads to an increased critical temperature of 
order T ~ ImK for the nuclear spins, which might be 
feasible in experiments. In Ref. [20| a C 13 carbon nan- 
otube was studied. By approximating the conduction 
electrons by a Luttinger liquid and treating the large ef- 
fective nuclear spins classically, the transition tempera- 
ture between a helically ordered (FM for finite systems) 
and unordered spin lattice was calculated.— It could be 
shown that a finite long-ranged Coulomb interaction is 
required to have a finite transition temperature^ which 
is consistent with the Mermin- Wagner Theorem 2 ^ and its 
recent extension.— Taking backaction effects of the nu- 
clear lattice on the electron spins into account increases 
the transition temperature by another order of magni- 
tude. This makes the KLM interesting for experiments, 
which are always performed at finite temperature. 

These developments motivate the study of the KLM 
in the presence of a finite Coulomb interaction between 
the itinerant electrons. The simplest extension to the 
KLM in terms of lattice models is the onsite Coulomb 
interaction U . In the case of half-filling a finite U leads 
to the opening of a spin and charge gap.— This work 
has been extended within a continuum Luttinger liquid 
approach to arbitrary fillings solved by bosonization^ 
Lattice effects have been accounted for by means of a 
phononic field and therefore there is no real lattice in- 
volved in those calculations. Still, the authors of Ref. 
find the interesting result of a shift of the phase boundary 
between FM and PM phase, as expected. 

In this paper, we use the density matrix renormaliza- 
tion group method^— (DMRG) to study ground state 
and dynamical properties of the one dimensional KLM 
for local spins with spin 1/2 including onsite and near- 
est neighbor Coulomb interaction. Our method benefits 



From ground state calculations we show that onsite 
Coulomb interaction lowers the value of J required for a 
transition from a PM to a FM ground state. For small 
n < 0.4 nearest neighbor Coulomb interaction V acts the 
same way on the magnetic order as U does. For n > 0.4 
they compete with each other. As a different sensor of 
magnetic order we utilize the static electron spin suscep- 
tibility. For the PM phase a peak at 2kp is expected 
(which diverges for L — > oo), while for the FM order 
a minimum at the smallest possible quasimomcntum q, 
which is finite for finite lattices, should emerge. This was 
stated similarly in Refs. I^l20l for small coupling constants 
J. 



Finally, we calculate the quasiparticle life-time broad- 
ening r + of an electron, its spin oriented in the opposite 
direction than that of all other electrons in the ground 
state. In Ref. [l?] it was shown in the FM phase and for 
electronic densities below half-filling that the effective in- 
teraction between spinpolaron states is weak proving that 
spinpolaron (spin-down) states are indeed well-defined 
quasiparticles with small life-time broadening T_ even in 
the presence of many electrons. However we will show 
here that the spin relaxation and decoherence rates will 
be dominated by the life-time broadening T + of the op- 
posite spin-up state, which is higher in energy. We will 
consider a single spin-up electron with quasimomentum k 
on top of the FM ground state of the Id KLM. Although 
this spin has the same direction as the underlying local 
spins and, thus, can not decay by direct exchange with 
the local spins, we find that T + is dominated by the ef- 
fective exchange interaction with the sea of spinpolaron 
spin-down states in the system. As a consequence, r + 
turns out to be much larger than T_ and dominates the 
spin relaxation as well as the spin decoherence rate (the 
pure dephasing term arising from the life-time broaden- 
ing r_ of the spin down spinpolaron state is negligible). 
We analyze the life-time broadening T + depending on J, 
U , n and the quasimomentum k and give explanations 
for the observations. Although the spin relaxation rate 
increases significantly in the presence of many electrons 
we will show in appropriate parameter regimes that the 
spin relaxation rate can be several order of magnitudes 
smaller in the FM phase compared to the PM phase. 



II. MODEL 

The Hamiltonian of the KLM with Coulomb interac- 
tion is sketched in Fig. Q] and defined as 
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FIG. 2. ( Color online) (a) Sketch of a configuration with three 
spinpolarons, each consisting of a delocalized spin singlet state 
with the local spins, (b) Sketch of a configuration with two 
spinpolarons and one spin up electron. 



where t is the hopping integral, L the lattice size, cfj the 
electron annihilation (creation) operator at site i with 
spin cr, J > the antiferromagnetic Kondo exchange 
coupling, S,; the spin operator of the local spin at site i, 
Si the spin operator of the conduction electron at site i, U 
the onsite Coulomb interaction constant, n ?(T = c\ a Ci a , V 
the nearest neighbor Coulomb interaction constant and 
Tii = + riii- All spins are considered to be spin 1/2. 
We define the filling by n — N/L, where N denotes the 
total number of itinerant electrons (n = 1 corresponds to 
half- filling). 



III. METHOD 
A. DMRG 

The DMRG method is a well established numerically 
exact method for the calculation of ground states, dy- 
namical properties and time evolution of one dimensional 
lattice systems. Our algorithm is formulated in a matrix- 
product language^! and makes use of Abelian, e.g. par- 
ticle number conservation (C/(l)) and non- Abelian, e.g., 
total spin conservation (SU(2)), symmetries. Depending 
on the symmetry sector, the use of SU (2) symmetries in 
addition to U(l) symmetries allows for computations up 
to 10 times faster. 



B. Ground states 

Calculating the ground state of a given system is syn- 
onymous to finding the symmetry sector with its corre- 
sponding quantum numbers, where the energy is mini- 
mal. The ground state phase diagram of the KLM is 
shown in Fig. [3] in dependence of the Kondo constant 
J and the filling n. Fixing J and n leaves the total 
spin quantum number S as the only free parameter, 
which distinguishes the order of the ground state, i. e., 
S = (L — N)/2 complies with FM order of local spins 
and S = with PM order. We choose SU{2) symmetry 
for the spin here, first for computational reasons and sec- 
ond it has the benefit that the states with different total 
spin quantum numbers arc non-degenerate in this case, 
whereas in U(l) symmetry a partial degeneracy in the 



total spin in the direction of quantization exists. Consid- 
ering Coulomb interaction in addition, we have another 
two variables that have to be fixed in advance and this 
means we have a quadruple of variables {n, J, U, V}, or 
a four dimensional phase diagram. 



C. Susceptibility 

We calculate the static electron spin susceptibility 
x(w = 0) by means of Green's functions and the applica- 
tion of dynamical DMRG^^ with GMRESM&. Details 
of our implementation can be found in Ref. Il7l . 

The definition of the spin susceptibility is 
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where H is the Hamiltonian given in Eq. (JTJ, |0) is the 
ground state of the system, and E the ground state en- 
ergy, rj is a finite artificial broadening factor, needed 
to avoid finite size effects^ and which can be choosen 
smaller with larger lattice size. 



D. Quasiparticle life-times 

In Ref. [I?] the quasiparticle life-time of the spinpolaron 
was calculated (cf. Fig. by evaluating the electronic 
Green's function in momentum and frequency space 



G ka (u> + irj) 



1 



u + irj — (eo(fc) — ix + T, a (k,u; + irj)) ' 



(3) 



where w is the energy, eo(fc) the free electron dispersion 
relation, fi the electrochemical potential (which does not 
play a role in the calculation of broadenings of spec- 
tral densities) and £<y(£;, u> + irj) the complex self-energy. 
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From the imaginary part of the self-energy, which is given 
by the broadening of the Lorentzian shaped peak in the 
spectral density A a (k,u) = — (l/V) ImGfc a (w) we can 
determine the quasiparticle life-time in dependence of 
all parameters. On the technical side, we use again the 
above mentioned GMRES method and calculate spectral 
densities as described in Ref. [T3- 

Basically, there exist four different scenarios for which 
the electronic quasiparticle life-time broadenings can be 
calculated assuming that in the FM ground state the 
local spins point up and the conduction electron spins 
point down (for large J the most dominant part of a spin- 
polaron state consists of a conduction electron pointing 
down with a small admixture of the spin up state plus a 
local magnon): 

• 1 In the FM phase for a spin down electron (cf. 
Fig. Hi); 

• 2 In the FM phase for a spin up electron (cf. 
Fig. [2b); 

• 3 and 4 are the corresponding cases for the PM 
phase. 

1 corresponds to the spinpolaron life-time broadening T_ 
and 2 to its natural counter part r+. 3 and 4 are iden- 
tical, since the spins in the PM ground state have no 
specific direction. 

In addition to Ref. [l?] we calculate here the life-time 
broadening T + . As shown in this paper this rate is very 
large in the presence of many electrons, T + 2> T_, and, 
as a consequence, dominates the spin relaxation and de- 
coherence rates, as can be understood from the following 
qualitative analysis. The two many-body spin states |±) 
depicted in Fig. [2] are not exact eigenstates but are ex- 
pected to be part of a sharp many-body continuum with 
long life-times. The spin down state |— ) is protected from 
magnon absorption and emission processes since the spin- 
polarons can lower their energy by the entanglement with 
the local spins in a singlet state. Only virtual processes 
and weak spinpolaron-spinpolaron interactions lead to a 
small broadening T_ of this state, as shown in detail 
in Ref. Ojl The spin- up state |+) is protected due to 
the spin polarization of the local spins. Due to effective 
exchange interaction between the spinpolarons and the 
spin-up electron mediated by the magnons, as discussed 
in detail in this paper in section lIV C[ this state has a life- 
time broadening ^> T_. Denoting the quasienergies 
of the two spin states by E± , we get a decay according 
to (±\e- lHt \±) ~ e -*E ± t e -(r ± /2)t_ Tq dcfine the gpin rc _ 

laxation and decoherence rates, we introduce pseudo-spin 
operators P z = (l/2)(|+)(+| - |-)(-|) and P± = \±)( T \. 
Using spin conservation, we obtain after a straigthfor- 
ward calculation that (P z (t)) = (l/2)\(+\ e - lHt \+)\ 2 , if 
the system is prepared at t = in the state |+), and 
(P+(t)) = {\/2){+\e- lHt \+)*{-\e- lHt \-) 1 if the system 
is prepared in the state (l/v / 2)(|— ) + |+)) intially. As 
a result we find for the two different initial preparations 
that (P z (t)) ~ e~ Tlt and (P+(t)) ~ e lAt e~ T2 \ where 



A = E + — E_ is the quasienergy splitting and the spin 
relaxation/decoherence rates are given by 

ri=r+ , r 2 = ir 1 + ir_ . (4) 

This result shows that the dominant part to T1/2 is given 
by the broadening T + of the spin- up state |+), whereas 
the broadening T_ of the spinpolaron state |— ) enters 
only into the pure dephasing term of longitudinal fluctu- 
ations and can be neglected. 

E. Dispersion relation 

The dispersion relation u a {k) can be constructed from 
the resonance of the single particle spectral density 
A a {k,ui) at uj = u a (k). The number of k values is re- 
stricted by the lattice size L. 

IV. RESULTS 

In nearly all cases we have choosen L = 48, which is 
suitable from two different points of view. First, physi- 
cally, we are especially interested in finite systems, which 
would more closely resemble, e.g., nanotubes in the real 
world. And second, from the point of view of compu- 
tational cost, it is not convenient to take larger systems 
into account, since we already needed up to 3000 DMRG 
states in some of the calculations, which is a large num- 
ber considering the number of executed calculations. All 
calculations are done with high computational precision, 
partly up to machine precision. We set t = 1 in all cal- 
culations. 



A. Phase diagram 

We will first investigate the influence of Coulomb in- 
teraction on the ground state of the Kondo lattice model. 
The phase diagram^ of the KLM (without Coulomb inter- 
action) is well established and shows two different phases, 
an FM and a PM one, see Fig. [3] The PM phase lies in 
the lower-right triangular of the phase diagram and for 
all other values of J and n < 1 the KLM has an FM 
ground state. Especially for N = 1 it was shown that 
the KLM is FM for any J ,12 As can be seen from Fig. [3l 
applying a finite onsitc Coulomb interaction shifts the 
phase boundary downwards for all values of n. This is 
consistent with the analysis of Ref. [2(| where a higher 
crossover temperature has been predicted in the pres- 
ence of Coulomb interaction. However, we note that the 
two mechanism are quite different. Whereas in Ref. 
the local nuclear spins have been treated quasiclassically 
due to their large effective spin, the present analysis is 
in the full quantummechanical regime of local spins with 
spin 1/2. Roughly speaking the present result is consis- 
tent with the Stoner picture of ferromagnctism, where 
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FIG. 3. (Color online) Phase boundary between the FM (up- 
per part) and the PM (lower part) ground state of the Kondo 
lattice model with L = 48 for three different cases of Coulomb 
interaction. 



FIG. 4. (Color online) Static electron spin susceptibility 
\(u) = 0) for a Kondo lattice model with L — 48, N = 6 and 
J = 0.15 with (dashed line) and without (solid line) Coulomb 
interaction. The thin vertical line marks in the PM phase. 



a finite Coulomb interaction leads to the preference of 
a fully spin-polarized state for the itinerant electrons. 
This state coincides with the qualitative picture of spin- 
polaron states pointing into the opposite direction of the 
local spins, see Fig. [2k, 

For finite nearest neighbor Coulomb interaction V we 
find the qualitatively different result, that the phase 
boundary is shifted downwards for n < 0.4 and upwards 
for n > 0.4 and therefore crosses the phase boundary of 
the KLM without Coulomb interaction. For small fill- 
ings this can be explained in the same way as for the 
onsite Coulomb interaction case. For filling n > 0.4 the 
electrons are relatively close to each other and therefore 
strongly influenced by V. The possibility to occupy the 
same site with two electrons of opposite spin does not 
lead to an increasing energy due to Coulomb interaction 
and increases the kinetic energy at the same time. There- 
fore, in this regime, the unordered state becomes more 
favorable. 

Summarizing, the onsite and nearest neighbor 
Coulomb interaction are concurring for small n < 0.4 
and behave competitively for large n > 0.4. These re- 
sults are pictured in Fig. [3j The solid blue line is the 
phase boundary of the non-interacting KLM. If Coulomb 
interaction is switched on, the phase boundary is lowered 
for all values of n (dashed dark blue line) . For U = and 
V finite, the phase boundary is lowered for small n and 
raised above the non-interacting case phase boundary for 
larger n. 
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FIG. 5. (Color online) Dispersion relation of a ^-electron in 
a KLM with J = 0.5, N = 4 and U = V = 0. 

spins for small J is 

Jrkky oc -X ± ( U! = 0, k, J, U). (5) 

Therefore the order of the local spin lattice should cor- 
respond to the absolute maximum of the static electron 
spin susceptibility. In Fig.Uwe show this for two extreme 
cases with L = 48 and N = 6. The first case (solid black 
line in the figure) with [7 = 0, V = has a PM ground 
state and shows the susceptibility in the non-interacting 
case. It has an absolute maximum at k = 2kp . This 
evidences that for the chosen set of parameters the state 
indeed orders paramagnctically in a RKKY like fashion. 
If Coulomb interaction is switched on with U = 10, V = 2 
(dashed brown line in the figure) the absolute maximum 
is at k = 0. In this case FM order becomes dominant. 



B. Susceptibilities 

For small J the order of the local spins manifests itself 
also in the static electron spin susceptibility. As was 
shown in Ref. the effective coupling between the local 



C. Dispersion relation 

We calculated the dispersion relation of a ^-electron in 
a KLM with L = 48, iV = 4, J = 0.5 and U = V = 0. 
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TABLE I. fc-dependence 


of relaxation rates for N = 


2 for different values of J and U. 
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k/(£4 


1)] 
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2 


3 
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.1= 


=0.5, U= 


=0 


0.00097 ± 0.00003 


0.00128 ± 0.00002 


0.00166 ± 0.00003 


0.00204 ± 0.00005 


.1= 


=0.5, U= 


=0.2 


0.00220 ± 0.00009 


0.00299 ± 0.00005 


0.00403 ± 0.00005 


0.0048 ± 0.0001 


.1= 


=0.3, U= 


=0 


0.00035 ± 0.00001 


0.000470 ± 0.000004 


0.00066 ± 0.00002 


0.00077 ± 0.00002 


J= 


=0.3, U= 


=0.2 


0.00146 ± 0.00006 


0.00198 ± 0.00004 


0.00280 ± 0.00007 


0.00331 ± 0.00008 
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FIG. 6. (Color online) Magnon density in the KLM with 
L = 4&, t = J = 1, U — V = in dependence of the quasimo- 
mentum. The number of electrons is varied between 2 and 8 
in steps of 2. 



The result is shown in Fig. [5] It shows a cosine shaped 
dispersion, which leads to the conclusion that the elec- 
tron behaves more or less like a free electron, only slightly 
affected by the presence of the local spin lattice. This can 
be explained by the fact that a f-electron can not flip its 
spin directly by an exchange process with a local spin due 
to spin conservation. In contrast, a .[-electron can do so, 
leading to the formation of spinpolarons, which can lower 
their energy by this process and obtain a larger effective 
mass leading to a sharper dispersion relation. However, 
as shown in the next section, the life-time broadening of 
f-electrons is generically larger than those of .[-electrons, 
since the decay processes for spinpolarons start in higher 
order in J than those for f-electrons. 



D. Quasiparticle life-times 



a) 

-L 

/ 

"/ 

/ 
"/ 
"/ 



b) 



</ = 

K2 — k\ 



1 

0.9 
0.8 
0.7 
0.6 
0.5 



j ftm "" 



It/4 n/2 




FIG. 7. (Color online) (a) Simplified itinerant electron band 
structure in fc-space. Light blue electrons on the left side are 
electrons initially in the ground state and electrons on the 
right side are additionally added to the ground state. The 
process shown correspond to a spin flip of the added electron 
at k = ki. After the spin flip, the electron has opposite spin 
with k — fa and has absorbed a magnon with q = fa — fa. 
(b) Accumulated magnon density pt ± as given in Eq. ([6]), for 
L = 48, N = 4 and J = 0.5. (c) Decay rates for L = 48, 
N — 4 and J = 0.5 in dependence of k. 



From the electronic spectral density A^(k, u>) we obtain 
the quasiparticle life-time broadenings T + in dependence 
of J, U, k and TV. As we calculate the Green's func- 
tion G-|-(fc, uj) in frequency space, we obtain two branches: 
The c|^- and the Cfc-j— branch, respectively. The first one 
corresponds to an additional electron placed in a certain 
/c-mode and interacting with the other electrons and the 
local spins. The second type addresses the spin up-part 
of the already existing electrons in the system. Therefore 



the two branches address two different sets of states in 
the spectrum of the Hamiltonian. Here we are interested 
in the first case only, since we would like to know, what 
happens to a spin up electron brought into the system in 
addition to the other electrons. 



7 



TABLE II. Relaxation rates in dependence of the electron number N, J and U. k is set to the lowest possible value k — n/(L+l). 
The given number of N is the number of electrons taken into account during the ground state calculations, i.e., the spin up 
electron is in addition to this number, (p) mark parameters, which correspond to the paramagnetic phase. 
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2 


3 


4 


6 


12 


.7 = 


=1.0, u=o 






0.00263 ± 0.00012 


0.00199 ± 0.00013 


0.00185 ± 0.00017 


0.00085 ± 0.00018 


J= 


=1.0, U=0.2 






0.00432 ± 0.00070 


0.00294 ± 0.00019 


0.00249 ± 0.00027 


0.0011 ± 0.0002 


.7 = 


=0.8, 17=0 






0.00188 ± 0.00008 


0.00130 ± 0.00009 


0.00087 ± 0.00009 


0.00213 ± 0.00024 


J= 


=0.8, 17=0.2 






0.00303 ± 0.00021 


0.00221 ± 0.00017 


0.00184 ± 0.00017 


0.00238 ± 0.00038 


J= 


=0.6, 17=0 






0.00115 ± 0.00005 


0.00078 ± 0.00005 


0.00081 ± 0.00007 


(p)0. 00508 ± 0.00031 


.7 = 


=0.6, U=0.2 






0.00239 ± 0.00012 


0.00156 ± 0.00012 






.1 = 


=0.5, U=0 


0.00104 ± 0.00002 


0.00097 ± 0.00003 


0.00082 ± 0.00003 


0.00062 ± 0.00004 


0.00066 ± 0.00005 




.7 = 


=0.5, U=0.2 


0.00233 ± 0.00004 


0.0022 ± 0.0001 


0.00205 ± 0.00008 


0.00126 ± 0.00011 


0.00141 ± 0.00013 




.1 = 


=0.5, U=0.4 


0.00142 ± 0.00021 


.7 = 


=0.5, U=0.6 


0.00245 ± 0.00018 


.7 = 


=0.5, U=0.8 


0.00384 ± 0.00045 


.7 = 


=0.3, U=0 


0.00041 ± 0.00001 


0.00035 ± 0.00001 


0.00033 ± 0.00001 


(p) 0.00144 ± 0.00010 


0.00150 ± 0.00020 




.7 = 


=0.3, U=0.2 


0.00158 ± 0.0003 


0.00146 ± 0.000016 


0.00123 ± 0.00008 


0.00084 ± 0.00007 






.7 = 


=0.1, U=0.0 


0.000030 ± 0.000001 




0.00004 ± 0.0000009 


(p) 0.00593 ± 0.00024 






.7 = 


=0.1, U=0.2 






0.00061 ± 0.00003 


(p) 0.00460 ± 0.00011 







1. Decay rate dependence on k 

In Tab. U we show decay rates of a spin up electron 
added to the N = 2 ground state. For all sets of U and J 
we find that the decay rate increases with increasing k as 
long as k is smaller than 2k p. Here we give an expla- 
nation considering momentum conservation and phase 
space arguments. In the FM ground state the lowest 
electronic orbitals in k space are occupied up to by 
the available electrons all with spin down. A state with 
wave vector k has quasimomcntum ±k due to the open 
boundary conditions. An additionally superimposed spin 
up electron with a certain wave vector k\ has to change 
to the state ky, > 2/cp in order to flip its spin, see Fig. [7^. 
This decay channel can only happen if a magnon is ab- 
sorbed with wave vector q = \ki ± k%\. Such magnons 
arc present in the ground state because each spinpolaron 
state consisting of a spin down electron with wave vec- 
tor k has a small admixture of spin up states with wave 
vector \k ± q\ and a local magnon in state q. Smaller 
values of k\ decreases the number of magnons with small 
wave vector q = \ki — k%\ to enable this process. This 
can be quantified by the magnon density per electron 
rn q = (S~Sq) /N (see Fig. and further by the accu- 
mulated magnon density 

= Y1 H m i ' ( 6 ) 

er=± q=\k 1 +ok 2 \ 

0<q<n , 2k F <k 2 <TT 

which is shown in Fig. [TJd and clearly states that the 
number of suitable magnons increases with increasing k\ 
even above 2kp until it falls off finally. This result qual- 
itatively reflects the decay rate for the spin up electron 
shown in Fig. Eh for a KLM with L = 48, N = 4, J = 0.5 
and U = V = 0. The decay rate first increases for small 
k as indicated by the accumulated magnon density. For 



values above 2/cp- the decay rate even surpasses the val- 
ues at 2k f until it decreases finally for larger values of 
k. We note that this is only a qualitative explanation 
since other decay channels involving absorption of many 
magnons are present as well. 

The discussed process for the decay of the spin up elec- 
tron is essentially an exchange process between a spin 
up electron in state k\ and a spinpolaron in state k. 
The spinpolaron provides the magnon with wave vec- 
tor q = \k\ ± to flip the spin up electron from state 
k\ — > k%, leaving the spinpolaron as a spin up electron 
in state |fc ± q\. As a result, by mediation of a local 
magnon, the spins of two electrons have been exchanged, 
whereas the local spin lattice is unaffected. This spin ex- 
change process is the essential process leading to a large 
life-time broadening of the spin-up electrons if many elec- 
trons are present in the system. In contrast, the spinpo- 
laron states have life-time broadenings, which are sev- 
eral orders of magnitude smaller compared to those of 
the spin-up states. The reason is that the spinpolaron- 
spinpolaron interaction is rather weak and can only be 
mediated via multi-magnon processes. 



2. Decay rate dependencies on U, J, N 

In this section we will explain how the quasiparticlc 
decay rate of the spin up state depends on U, J and N 
and why the found tendencies are to be expected. The 
results for these cases are shown in Tab. |TTJ 

Let us first consider the J-dependency. Picking one of 
the columns and considering only one of the two U -values 
we immediately recognize that the decay rate shrinks 
with decreasing J. The exchange strength J determines 
the time scale on which spins will flip, therefore with 
decreasing J flipping will be suppressed and the rate de- 
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creases. We note that this is different for the decay rate of 
the spinpolaron, where an increasing J stabilizes each po- 
laron and makes it insensitive to interactions with other 
electrons. For small J close to or even in the PM phase 
the decay rate of the spin up state increases notably, see 
A = 4. This is natural, since in a paramagnetically or- 
dered system many additional decay channels will open 
up. 

Considering the U dependence we find that with in- 
creasing U the rate increases in most cases. In section 
IIV Al we have found that an onsite Coulomb interaction 
has the tendency to order the local spins ferromagncti- 
cally. The additional spin up electron tries to align paral- 
lel to the other electrons to minimize interaction energy 
from the Coulomb potential. This infers a larger decay 
rate, if U becomes larger. Therefore this tendency here 
complies with the influence of the onsite Coulomb inter- 
action found above. Only when a finite U triggers the 
crossover from the PM to the FM phase, the rate de- 
creases with increasing U, sec A = 4 and J = 0.3. This 
is obvious since in the PM phase the phase space argu- 
ments presented in section HV D H are no longer valid and 
many more decay channels are possible. 

If we increase the number of electrons N in the system 
and keep the quasimomentum k fixed we find that the 
rates decrease with increasing N, for small N deep in 
the FM phase. This can be explained analog to the dis- 
cussion in section llVD 11 In the ground state, all initially 
available electrons fill the spinpolaron-band successively 
up to 2kp mainly in the spin down state. An additional 
spin up electron can be added to any fc-mode. In Tab. HE 
we considered the lowest state k = tt/(L + 1) in all cases. 
Considering one of the rows the electron number is in- 
creased from left to right and with each electron more in 
the ground state the respectively next higher fc-mode is 
occupied by this additional electron. As a consequence, 
as shown in section llVD 11 by increasing A we decrease 
the number of magnons suitable for scattering processes 
and therefore the decay rate has to decrease. However, 
in competition to this effect, increasing A" means also 
approaching the PM phase. Then we expect that differ- 
ent and also more decay channels open up, which should 
lead to an increasing decay rate. This can be seen in 
Tab. HE for J = 0.8 between A = 6 and A = 12. We 
have also calculated lifetimes for A = 7,9,10,11 (not 
shown), showing that the decay rates are monotonically 
increasing with increasing N for large N. For values of 
N close to half-filling of the conduction band and large 
values of J, such that we can switch between PM and 
FM phase, we find decay rates of the order of 0.01. As 
a consequence, the decay rate depends nonmonotonically 
on A, it decreases for small values of A deep in the FM 
phase and increases for larger values of N when the PM 
phase is approached. 

Nonetheless we find the sweet spot of the system by 
decreasing the number of electrons going from A = 4 to 
A = 3 electrons at J = 0.1. There we find that the de- 
cay rate of the spin up electron decreases by two orders of 



magnitude when comparing the rates in the PM and FM 
phase. Still it is important to note that a minimum num- 
ber of electrons in the system is important to maintain 
the FM order, especially at finite temperatures. 



V. DISCUSSION 

In this work we discussed the phase diagram and the 
spin relaxation properties of the Id spin-1/2 Kondo lat- 
tice model with Coulomb interaction. We found that a 
finite onsite or nearest neighbor interaction favors a FM 
order of the local spin lattice for small enough electronic 
densities. This gives further strong support to the anal- 
ysis of Refs. I3l20l where similiar results have been found 
in 2d semiconductor systems and C 13 carbon nanotubes. 
It provides a pathway to achieve a spontaneous and full 
polarization of the nuclear spins by lowering the temper- 
ature below the critical one. This configuration is desir- 
able for applications in quantum information processing, 
since it reduces the spin relaxation and decoherence rates 
of the electronic spins. It is important to notice that a 
finite crossover temperature can only be expected, if the 
density of electrons is finite. Thus, many electrons are 
necessary to achieve the FM state. Once the FM state is 
achieved, one can in principle perform quantum informa- 
tion processing by realizing quantum dots with external 
gates on time scales which are small compared to the time 
the nuclear spins need to return to the PM phase. If this 
is possible one can effectively realize a system consisting 
of one single electron AT = 1 in contact with a ferromag- 
netically ordered nuclear spin lattice. In this case the spin 
up state and the spinpolaron are exact eigenstatcs, i.e., 
the ideal situation with T± = is achieved. In this paper 
we discussed the spin relaxation properties for N > 1, i.e. 
we analysed the question whether the spins in a many- 
body system could possibly be used as candidates for 
spin qu-bits. In Ref. [l?] we already found that spinpo- 
larons arc indeed very long living states, indicating that 
the spinpolaron-spinpolaron interaction is rather weak. 
However, in this paper we found that the spin up state is 
strongly influenced by exchange interaction between the 
spin up and spinpolaron states. This exchange process 
does not require any finite energy and, therefore, can not 
even be suppressed by application of a finite magnetic 
field. We analysed in detail the dependence of r + on the 
Coulomb interaction U, the exchange interaction J, the 
particle number A" and the quasimomentum k. In the 
FM phase we found that the rate decreases for smaller 
values of U, J, k, and larger values for N, unless we ap- 
proach the PM phase. For appropriate parameter sets 
we have shown that the life-time of spin up states can 
be two orders of magnitude larger in the FM phase than 
in the PM phase. However, compared to the life-time of 
spin down spinpolaron states, their life-time is orders of 
magnitudes smaller, regardless of the chosen parameter 
regime in the FM phase. 
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